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Fermi surface topology of the two-dimensional Kondo lattice model: a dynamical 

cluster approximation approach 

L. C. Martin, M. Bercx and F. F. Assaad 

Institut filr Theoretische Physik und Astrophysik, 

Universitdt Wiirzhurg, Am Hubland, D- 97074 Wilrzhurg, Germany 

We report the results of extensive dynamical cluster approximation calculations, based on a 
quantum Monte Carlo solver, for the two-dimensional Kondo lattice model. Our particular cluster 
implementation renders possible the simulation of spontaneous antiferromagnetic symmetry break- 
ing. By explicitly computing the single-particle spectral function both in the paramagnetic and 
antiferromagnetic phases, we follow the evolution of the Fermi surface across this magnetic tran- 
sition. The results, computed for clusters up to 16 orbit als, show clear evidence for the existence 
of three distinct Fermi surface topologies. The transition from the paramagnetic metallic phase to 
the antiferromagnetic metal is continuous; Kondo screening does not break down and we observe 
a back- folding of the paramagnetic heavy fermion band. Within the antiferromagnetic phase and 
when the ordered moment becomes large the Fermi surface evolves to one which is adiabatically 
connected to a Fermi surface where the local moments are frozen in an antiferromagnetic order. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27.+a, 73.43.Nq, 75.30.Kz, 87.15.ak 



INTRODUCTION 



Heavy-fermion systems^ are characterized by a hier- 
archy of distinctive energy scales^^^. The Kondo scale, 
Tk oc e"^/'^ with bandwidth W and superexchange J, 
marks the screening of local magnetic moments. This 
screening is a many-body effect which entangles the spins 
of the conduction electrons and local moments^. Below 
the coherence temperature, which is believed to track the 
Kondo scale^^, the paramagnetic (PM) heavy-fermion 
liquicP emerges and corresponds to a coherent, Bloch 
like, superposition of the screening clouds of the individ- 
ual magnetic moments. Even in the Kondo limit, where 
charge fluctuations of the impurity spins are completely 
suppressed, this paramagnetic state is characterized by a 
large Fermi surface with Luttinger volume including both 
the magnetic moments and conduction electron d^^^^^ l 
The coherence temperature of this metallic state is small 
or, equivalent ly, the effective mass large. 
Kondo screening competes with the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction, which indirectly 
couples the local moments via the magnetic polarization 
of the conduction electrons. The RKKY energy scale is 
set by J'^Xc{'^i^ = 0) where Xc corresponds to the spin 
susceptibility of the conduction electrons^. 
The competition between Kondo screening - favoring 
paramagnetic ground states - and the RKKY interaction 
- favoring magnetically ordered states - is at the heart of 
quantum phase transitions^^ ^^, the detailed understand- 
ing of which is still under debate (for recent reviews see 

Ref. HMEl). 

Here, two radically different scenarios have been put for- 
ward to describe this quant um ph ase transition. In the 
standard Hertz-Millis picture^^^^, the quasi-particles of 
the heavy-fermion liquid remain intact across the tran- 
sition and undergo a spin-density wave transition. In 
particular, neutron scattering experiments of the heavy- 



fermion system Cei_a;Laa^Ru2Si2 show that fluctuations 
of the antiferromagnetic order parameter are responsi- 
ble for the magnetic phase transition and that the tran- 
sition is well understood in terms of the Hertz-Millis 
approaclP^. 

On the other hand, since many experimental observa- 
tions such as the almost wave vector independent spin 
susceptibility in CeCue-ajAuajl^, or the jump in the low- 
temperature Hall coefficient in YbRh2Si2^^ are not ac- 
counted for by this theory alternative scenarios have been 
put forwarcP^H^. In those scenarios, the quantum crit- 
ical point is linked to the very breakdown of the quasi- 
particle of the heavy-fermion state^^ ^^, and a topological 
reorganiza tion o f the Fermi surface across the transition 
is expected^^ ^^. 

Recent experiments on CelnJ^ or CeRhi_a^Coa^InJ^ 
show that a change in Fermi surface (FS) topology must 
not necessarily occur only at the magnetic order-disorder 
quantum critical point (QCP). In fact, even in YbRh2Si2 
it has since been shown that the Fermi surface recon- 
struction can be shifted to either side of the QCP via 
application of positive or negative chemical pressure^. 

In this paper, we address the above questions through 
an explicit calculation of the Fermi surface topology in 
the framework of the Kondo lattice model (KLM). In its 
simplest form the KLM describes an array of localized 
magnetic moments of spin 1/2, arising from atomic /- 
orbitals, that are coupled antiferromagnetically (AF) via 
the exchange interaction J to a metallic host of mobile 
conduction electrons. 

We present detailed dynamical cluster approximation 
(DCA) calculations aimed at the investigation of the 
KLM ground state. For the simulations within the mag- 
netically ordered phase, we have extended the DCA to al- 
low for symmetry breaking antiferromagnetic order. We 
map out the magnetic phase diagram as a function of 
J/t and conduction electron density ric^ with particular 
interest in the single-particle spectral function and the 



evolution of the Fermi surface. The out hue is as fohows. 
The model and the DCA implementation is discussed in 
Sec. [TTj Results for the case of half-band filling and hole- 
doping are discussed in Sec. Ill and IV Section [V| is 



devoted to a summary. This paper is an extension to our 
previous work, where part of the results have already 
been published^. 



II. MODEL HAMILTONIAN AND DYNAMICAL 
CLUSTER APPROXIMATION 

The Kondo lattice model (KLM) we consider reads 

iJ = ^(e(k)-M)4,.Ck,a + jESi-Sf ■ (1) 

k,cr i 

The operator cj^ ^ denotes creation of an electron in 
a Bloch state with wave vector k and a z-component 
of spin (J =t, i- The spin 1/2 degrees of freedom, 
coupled via J > 0, are represented with the aid of 
the Pauli spin matrices a by S? = | J^s s' ^i s^s^^'^i^s' 
and the equivalent definition for Sj using the localized 
orbital creation operators /j ^. The chemical potential 
is denoted by /i. The definition of the KLM excludes 
charge fiuctuations on the /-orbitals and as such a strict 
constraint of one electron per localized /-orbital has to 
be included. For an extensive review of this model we 
refer the reader to Ref. |36l 

Particle-hole symmetry at half-filling is given if hopping 
is restricted to nearest neighbors on the square lattice and 
the chemical potential is set to zero. We introduce a next- 
nearest neighbor hopping with matrix element t' to give 
a modified dispersion e(k) = —2t [co^ikx) + co^{ky)] — 
2t' [cosikx + ky) -h co^{kx — ky)]. As we will see, the 
suppression of particle-hole symmetry due to a finite 
value of t' leads to dramatic changes in the scaling of the 
quasi-particle gap at low values of J /t and at half-band 
filling. We have considered the value t' /t = —0.3. This 
choice guarantees that the spin susceptibility of the host 
metallic state at small dopings away from half-filling is 
peaked at wave vector Q = (7r,7r). Thus, antiferromag- 
netic order as opposed to an incommensurate spin state 
is favored. 



To solve the model we have used the DCA^^ ^^ ap- 
proach which retains temporal fiuctuations and hence 
accounts for the Kondo effect but neglects spatial fiuc- 
tuations on a length scale larger than the cluster size. 
The approach relies on the coarse-graining of momentum 
space, and momentum conservation holds only between 
the k-space patches. By gradually defining smaller sized 
patches the DCA allows to restore the k-dependency of 
the self-energy. 

The standard formulation of DCA is naturally transla- 
tionally invariant. To measure orders beyond translation 
invariance the DCA can be generalized to lattices with a 
supercell containing A/'^-unit cells of the original lattice^^. 



A general unit cell is addressed by x = x + r^, x denoting 
the supercell and r^ the relative points with ji = l...Nu- 
In this work we opted for N^ = 2 which is the minimum 
requirement for the formation of antiferromagnetic or- 
der. The reduced Brillouin zone (RBZ) is then spanned 
by bi = 7r(l, 1) and b2 = 7r(l, —1) and we have set the 
lattice constant to unity. This supercell as well as the 
RBZ (also coined magnetic Brillouin zone (MBZ)) are 
plotted in Fig. [l] 

The self-energy and Green function are to be understood 
as spin dependent matrix functions, with an index for the 
unit cell within the supercell as well as an orbital index 
for the c- and /-orbitals in each unit cell. The k-space 
discretization into Np patches with momentum conserva- 
tion only between patches yields the coarse-grained lat- 
tice Green function 



k 



(2) 



and ensures that the self-energy is only dependent on the 
coarse-grained momentum K: S = I]cr(K,ia;). Here, the 
reciprocal vector K denotes the center of a patch and 
the original k vectors are given by k = K + k. The 
self-energy is extracted from a real-space cluster calcu- 
lation with periodic boundaries yielding the quantized 
K values. Let ^o,cr(K,icj) be the bath (non-interacting) 
Green function of the cluster problem and ^cr(K, iuj) the 
full cluster Green function. Hence, S = Qq — Q~^ and 
self-consistency requires that: 



aa(K,ia;) = (ao"i(K,ic^) - S,(K,ic^))- 



(3) 



N 



f^(Go-i(K + k,ia;)-S.(K,ic^))-^ 



The non-interacting lattice Green function is denoted by 
Go,a(K + k,ia;). 

To summarize, our implementation of DCA approximates 
the Fourier space of a lattice with two-point basis by 
patching. Regarding the KLM this leads to real-space 
clusters of size Np , each encompassing 2 x Np c- and /- 
orbitals. In the present work cluster sizes Np = 1 and 
Np — 4 are considered (F: 



Fig. [T]). 

ion G^(k,icjn) 



The lattice green function Gcr(k, icj^) with k G MBZ 
can equally be expressed as ^cr(k, icj^) in an extended 
Brillouin zone scheme, with k G BZ and k = k + vf\h\ -|- 



1 



A^«=2 



^.(k,ia;.) = - ^ e^^(r-r^)[G.(k,ic..)]^, . (4) 

In order to implement the DCA one has to be able to 
solve, for a given bath, the Kondo model on the effec- 
tive cluster. The method we have chosen is the auxil- 
iary field Quantum Monte Carlo (QMC) version of the 
Hirsch-Fye algorithm, following precisely the same real- 
ization as in Ref. [39l The implementation details of the 
self-consistency cycle can be found in Ref. |40l The per- 
formance of the QMC cluster solver has been enhanced 
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FIG. 1. (Color online) Definition of the antiferromagnetic unit 
cell (a) and DCA patching of the Brillouin zone (b). The lat- 
tice vectors ai = (1,1) and 3.2 = (1,-1) connect the AF unit 
cells. The two inequivalent c-orbitals are denoted by filled 
(empty) circles and the localized /-orbitals with arrows. The 
reciprocal-space lattice vectors bi and b2 span the magnetic 
Brillouin zone. The MBZ is uniformly discretized in Np — 4 
patches, corresponding to real- space clusters each containing 
Np magnetic unit cells. The color coding refers to constant 
momentum dependency of the self-energy and cluster Green 
function. The dashed square denotes the extended Brillouin 
zone. 



by almost an order of magnitude by implementing the 
method of delayed updates"^^. During the QMC Markov 
process the full updates of the Green function are delayed 
until a sufficiently large number of local changes in the 
auxiliary field have been accumulated. This enhances the 
performance of the whole QMC algorithm since consid- 
erable CPU time is spent in the updating section. To ef- 
ficiently extract spectral information from the imaginary 
time discrete QMC data a stochastic analytic continua- 
tion scheme is employed^. 

Ground state properties of the model are accessed by ex- 
trapolation of the inverse temperatures (3 to infinity. This 
is a demanding task since all relevant energy scales, the 
Kondo scale, the coherence scale and the RKKY scale, 
become dramatically smaller with decreasing J/t. In the 
paramagnetic phase and below the coherence scale a clear 
hybridization gap is apparent in the single particle spec- 
tral function and we use this criterion to estimate the 
coherence temperature. Since the computational time 
required by the QMC cluster solvei^^ increases propor- 
tional to (PNuNp)^^ this limits us in the resolution of 
energy scales to values of J/t > 0.8. It is important 
to note that for small dopings away from half-band fill- 
ing and for the considered cluster sizes, the negative sign 
problem is not severe and hence is not the limiting factor. 



III. THE HALF-FILLED KLM 

We consider the KLM at half- filling in two cases: either 
with or without particle- hole symmetry. That is t' /t = 
and t' /t = —0.3 respectively. At t' = lattice QMC sim- 



ulations do not suffer from the negative sign problem, 
and it is well established that the Kondo screened phase 
gives way to an AF ordered phase at Jc/t = i.4 J39 | 44 |^ 
This magnetic order-disorder transition occurs between 
insulating states and it is reasonable to assume that it be- 
longs to the three-dimensional 0(3) universality clasa^. 
Here, the dynamical exponent takes the value of unity 
such that the correlation lengths in imaginary time and 
in real space are locked in together and diverge at the 
critical point. Due to the very small cluster sizes consid- 
ered in the DCA it is clear that we will not capture the 
physics of this transition. In fact as soon as the correla- 
tion length exceeds the size of the DCA cluster, symme- 
try breaking signaled by a finite value of the staggered 
moment 



mi 



2Np^ 



(n/. 



^i,; 



)e- 



iQi 



(5) 



sets in and mean-field exponents are expected. The nor- 
malization is chosen such that the staggered magnetiza- 
tion of the fully polarized state takes the value of unity. 
The above quantity is plotted in Fig. |2] at t' /t = and 
Np = 1 As apparent, the critical value of J/t overes- 
timates ( Jc/t ^ 2.1.) the exact lattice QMC result. 
Switching off nesting by including a finite value of t' /t 
shifts the magnetic order-disorder transition to lower val- 
ues of J/t. 

To improve on this result, we can systematically enhance 
the cluster size. However, our major interest lies in 
the single-particle spectral function. As we will see be- 
low, and well within the magnetically ordered phase, this 
quantity compares very well with the exact lattice QMC 
results. 



A. Single-particle spectrum - t' /t = 

We have calculated the single-particle spectral func- 
tion A^^(k,a;) = —^ X^cr-"-^ ida^O^^^)] of ^^^ conduction 
electrons along a path of high symmetry in the extended 
Brillouin zone. In all our spectral plots, on the energy 
axis, CO values are given relative to the chemical potential 
/i. 

In Fig. [3] we plot the spectrum for J/t = 2.2 which was 
seen to be in the paramagnetic phase (see Fig. [2|. The 
lower band runs very flat around k = (7r,7r), with rel- 
atively low spectral weight (note the logarithmic scale 
of the color chart) in comparison to the other parts of 
the band which are mostly unchanged from the non- 
interacting case. This feature is associated with Kondo 
screening of the impurity spins and the resultant large 
effective mass of the composite quasi-particles. Since no 
band crosses the Fermi energy {uj/t = 0) we classify this 
region of parameter space as a Kondo insulator. The ob- 
served dispersion relation is well described already in the 
framework of the large- A/" mean-field theory of the Kondo 
lattice modeF. At the particle- hole symmetric point this 
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FIG. 2. The staggered magnetization m{ of the local mo- 
ment spins in the ground state as a function of coupling J/t 
with next-nearest neighbor hopping t'/t = —0.3 and t'/t = 0, 
respectively. The low temperature limit is reached by per- 
forming simulations at various inverse temperatures /3 . Con- 
vergence to the ground state has been achieved for the fol- 
lowing inverse temperatures: /3t = 20 (2.5 > J/t > 1.7), 
/3t = 40 (1.6 > J/t > 1.1), I3t = 80 {J/t = 1) and /3t = 100 
(J/t = 0.8). 



approximation, which in contrast to the DCA+QMC re- 
sults of Fig. [3] neglects the constraint of no double occu- 
pancy of the /-orbitals, yields the dispersion relation of 
hybridized bands: 



E±(k) = - e(k)±Ve(k)2 + A2 



(6) 



Within the mean-field theory the quasi-particle gap (see 
Eq. m), Aqp (X A^, tracks the Kondo scale. 

At J/t = 2.0 we have measured a non- vanishing stag- 
gered magnetization m( = 0.244 ± 0.001. The spectral 
function (Fig. [4| for this lightly AF ordered simulation 
now includes additional low-energy band structures: In 
the upper band around k = (tt, tt) and in the lower band 
around k = (0,0). These shadow bands arise due to the 
scattering of the heavy quasi-particle off the magnetic 
fluctuations centered at wave vector Q = (7r,7r). 

At J/t = 1.2 we make a comparison with the spectrum 
obtained via QMC lattice simulations using the projec- 
tive auxiliary field algorithm of Ref. [39] for the same 
parameters. We see in Fig.[5]the QMC lattice simulation 
results^^ on the left and our DCA result on the right, 
both only for the photoemission spectrum {u/t < 0). The 
agreement of the DCA spectrum with the Blankenbecler- 
Scalapino-Sugar (BSS) result is excellent. As in the BSS 
lattice calculations and for the particle-hole symmetric 
case, the conduction (valence) band maximum (mini- 
mum) is located at k = (tt, tt) ( k = (0, 0)) for all consid- 
ered values of J/t. 

It is noteworthy that already with the smallest possible 
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FIG. 3. (Color online) Single-particle spectral function 
A^^(k, cj) at half-filling (ric = 1) and with particle-hole sym- 
metry (t'/t = 0) close to the magnetic phase transition on the 
paramagnetic side, J/t = 2.2. 
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FIG. 4. (Color online) Single-particle spectral function 
A^^(k,c<;) at half- filling {ric = 1) and with particle- hole sym- 
metry {t'/t = 0) close to the magnetic phase transition on the 
AF side, J/t — 2.0 (a) and closeup (b). The onset of magnetic 
order and concomitant emergence of the MBZ is signaled by 
the appearance of shadow bands. The staggered magnetiza- 
tion for this point is measured to be m{ = 0.244 di 0.001. 
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FIG. 5. (Color online) Comparison of the DCA single-particle 
spectral function A^^(k,a;) with the T = spectral function. 
Data (a) taken from Ref. |39]which used a projective auxiliary 
field Monte Carlo (BSS) method for the KLM on a 12 x 12 
lattice and our DCA results (b) at /3t = 40, both for the same 
parameter set - t'/t = 0, ric = 1, J/t = 1.2. 



cluster capable of capturing AF order, A/p = 1, we are 
able to produce a single-particle spectrum which is essen- 
tially the same as the lattice QMC result. This confirms 
that the DCA is indeed a well suited approximation for 
use with the KLM: the essence of the competition be- 
tween RKKY-mediated spacial magnetic order and the 
time-displaced correlations responsible for Kondo screen- 
ing is successfully distilled to a small cluster dynamically 
embedded in the mean-field of the remaining bath elec- 
trons. 



B. Single-particle spectrum - t' /t 



-0.3 



We now set t' /t = —0.3, but remain at half-filling via 
careful adjustment of the chemical potential fi. With de- 
creasing J /t we can divide the phase diagram into four 
regions on the basis of the characteristic spectrum in each 
region. For J/t = 2.5 (Fig.JGja)) the spectrum is qualita- 
tively identical to the PM phase in the t' /t = case such 
that the model is a Kondo insulator with an indirect gap. 
The minimum of the valence band lies at k = (0, 0) and 
the maximum of the conduction band is at k = (tt, tt). 
For J/t = 1.8 (Fig. [OJ^b)) the system is weakly AF or- 
dered, and the spectrum refiects this in the develop- 
ment of shadow bands, but is otherwise qualitatively the 
same as before. At J/t = 1.4 (Fig. [6]^c)), well inside 
the AF phase, the shadow bands are more pronounced 
and now the minimum of the valence band has shifted to 
k = (0,7r). The final plot in the series (Fig. JGJ^d)), with 
J/t = 1.0 shows that the form of the lower band has also 
changed such that the maximum of this band now lies at 
k=(7r/2,7r/2). 

In Figs. IrFa-h) we show the evolution of this lower band 
between k = (0, 0) and k = (tt, tt) with an enlargement 
of the energy axis around the Fermi energy. The local 
minimum energy dip at k = (7r/2, 7r/2) becomes less pro- 
nounced as the heavy-fermion band flattens until by the 
time we reach J/t = 1.2 this dip has turned into a bump 
such that the lower band maximum has shifted from 
k = (7r,7r) to k = (7r/2,7r/2). This important change 
in band structure appears to be continuous. For smaller 
coupling, J/t = 1.1, 1.0, and 0.8 the lower band maxi- 
mum remains at k = (7r/2,7r/2) and becomes more pro- 
nounced. Assuming a rigid band picture, this evolution 
of the band structure maps onto a topology change of 
the Fermi surface at small dopings away from half filling. 
We will see by explicit calculations at finite dopings that 
this topology change indeed occurs. 
The position of the minima and maxima in the valence 
and conduction bands of the spectra for t' /t = —0.3 are 
summarized schematically in Fig. |8] 



C. The quasi-particle gap 

The quasi-particle gap corresponds to the energy dif- 
ference between the top of the conduction band and the 



bottom of the valence band. To be more precise, 

2A,p = min «+i(k) - ^o^) - max {E^ - E^-\k)) . 

(7) 
Here, Eq ^(k) corresponds to the ground state energy 
of the N ±1 particle system with total momentum k and 
£^0^ = mink£^o^(k). Technically, the energy differences 
are at best extracted from the low temperature single- 
particle Green functions: 

lim (c(T)k4) oc e-«+'(k)-<-^) ,r » 1 

/3^-oo 

lim (4(T)ck> oc e-(<-<-'(k)-^) ,r « -1 . 

In the large J/t limit, each impurity spin traps a con- 
duction electron in a Kondo singlet. The wave function 
corresponds to a direct product of such Kondo singlets 
and the quasi-particle gap is set by the energy scale 4 J/3 
required to break a Kondo singlet. Hence, in the large J/t 
limit the quasi-particle gap does not depend on the de- 
tails of the band structure. The situation is more subtle 
in the weak-coupling limit, where the underlying nest- 
ing properties of the Fermi surface play a crucial role. 
Assuming static impurity spins locked into an antiferro- 
magnetic order, (Sj ) = ^e^e*^'^ one obtains the single 
particle dispersion relation 



^±(k) = e+(k)±Je2_(k) + (J/4)2, 



(8) 



with e±(k) = ^^ ^ ^^ — —. In one dimension and at 
tyt = nesting leads to e+(k) = and magnetic order 
opens a quasi-particle gap set by Aqp = J/4. Both the 
DCA results presented in Fig. [9]as well as the BSS results 
of Ref. 39 support this point of view. In one-dimension, 
this scaling of the quasi-particle gap is also observed^. 
At t' /t = —0.3 nesting is not present and magnetic order- 
ing can only partially gap the Fermi surface. Since the 
DCA results of Fig. ([9| support a finite quasi-particle 
gap in the magnetically ordered phase, the above frozen- 
spin Ansatz fails and points to one of the main result of 
our work, namely that within the magnetically ordered 
phase Kondo screening and the heavy quasi-particles are 
still present. Even though we cannot track the functional 
form of the gap at t' /t = —0.3 and at weak couplings, we 
interpret the DCA results in terms of a quasi-particle 
gap which tracks the Kondo scale from weak to strong 
coupling as in the large- A^ mean-field calculations^. 



IV. THE HOLE-DOPED KLM 

In this section we concentrate on the hole-doped KLM 
at t' /t = —0.3. We will first map out the magnetic phase 
diagram in the doping versus coupling plane and then 
study the evolution of single-particle spectral function 
from the previously discussed half-filled case to the heav- 
ily doped paramagnetic heavy-fermion metallic state. 
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FIG. 6. (Color online) Single-particle spectral functions A^^(k, cj) at half-filling {ric = 1) with t' /t = —0.3. As the Kondo cou- 
pling J/t is decreased, the spectrum changes from a paramagnetic Kondo insulator with indirect gap (a) to an antiferromagnetic 
ordered state (d). 



A. Magnetic phase diagram 



Fig. 10 shows the staggered magnetization as a func- 
tion of conduction electron density for J/t = 0.8, 1.0, 1.2 
and 1.4. The magnetically ordered state found at half- 
filling initially survives when doping with holes. At all 
coupling values a continuous magnetic phase transition is 
observed with the AF order decreasing gradually as the 
system is doped and vanishing smoothly at a quantum 
critical point. The results are summarized in the mag- 
netic phase diagram shown in Fig. 11 With decreasing 



values of J/t and increasing dopings the RKKY interac- 
tion progressively dominates over the Kondo scale and 
the magnetic metallic state is stabilized. Checks were 
made by varying the temperature, 1//3, of the simula- 



tions to ensure that the results can be considered to be 
ground state. Below J/t = 0.8 we found the coherence 
scale to be too low to guarantee convergence. As pointed 
out previously, the limiting computational factor is the 
cubic scaling of the Hirsch-Fye^ algorithm, {(3Np)^ and 
not the negative sign problem. 



B. Single-particle spectrum and topology of the 
Fermi surface 



We have calculated and followed the evolution of the 
single-particle spectrum, plotted in an extended Bril- 
louin zone scheme. To begin with we set J/t = 1 and 



show results in Fig. 12 for the spectral function, as it 
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FIG. 7. (Color online) Closeup for low energies around the Fermi energy (given by oj/t = 0). Development of the single-particle 
spectral function A^^(k, cj) at half-filling with t^t — —0.3 which shows a shift of the maximum in the lower band from k = (tt, tt) 
to k = (7r/2,7r/2) with decreasing J/t. 
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FIG. 8. (Color online) Schematic representation of the posi- 
tion of the upper band minima (top row) and lower band max- 
ima (bottom row) of the conduction electron single-particle 
spectral function with next-nearest neighbor hopping t' /t — 
—0.3 and for half-filling. The squares represent the first Bril- 
louin zone (for the PM phase) or the extended Brillouin zone 
(for the AF phase) , with the bottom left corner and top right 
corner of each square given by k = (— tt, — tt) and k = (tt, tt), 
respectively. 
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FIG. 9. The quasi-particle gap Aqp/t as a function of Kondo 
coupling J/t. At t' = — 0.3t, the magnetic order-disorder 
transition takes place at Jc/t ^ 1.85. In the coupling range 
0.8 < J/t < Jc/t Kondo screening coexists with magnetism 
and is at the origin of the quasiparticle gap. 



is here that a change in the Fermi surface topology be- 
comes evident. Starting from the previously discussed 
half-filled case (Fig. 12 ^a)) a rigid band approximation 
produces hole pockets around the k = (±7r/2, ±7r/2) 
points in the Brillouin zone. This is confirmed by ex- 
plicit calculations at small dopings away from half-filling 
as demonstrated in Fig. 12 ^b). With further doping this 
low energy band flattens out progressively becoming al- 
most flat by Tie = 0.926 (Fig. fl2|d)). Doping further 



gives rise to a Fermi surface with holes centered around 
k = (±7r,±7r). Since at Uc = 0.908 (Fig. ^e)) and 
Tie = 0.898 (Fig. Il2|f)) we still have non-zero magne- 
tizations of m( = 0.340 and 7ti{ = 0.245, respectively, 
shadow features centered around k = (0,0) are expected. 
The weight of those shadow features progressively di- 
minishes as the staggered magnetization vanishes. In 
particular, in Fig. fl2[g) the magnetization is very small 
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FIG. 10. The staggered magnetization m{ of the /-electrons 
as a function of Uc at different constant couphngs J/t. 
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FIG. 11. (Color online) Ground state magnetic phase dia- 
gram of the hole-doped KLM showing simulation results for 
the staggered magnetization ml (color-coded) as a function 
of coupling J/t and conduction electron occupancy Uc. Tri- 
angles: PM region, large FS. Squares: AF, large FS. Circles: 
AF, small FS. Here t' /t — —0.3 and the calculations are car- 
ried out with the A/p = 1 cluster. Below, the FS topologies 
corresponding to the numbered regions are shown schemati- 
cally. 



{ml = 0.123) and indistinguishable from zero in the final 
figure of the series. 

Across the transition from Uc = 1 to nc = 
0.856 the quasiparticle weight Z(k = (7r,7r)) = 
|(^^+-^|c|^_. . ^l^o^)P as obtained from the behaviour 
of ^^^(k = (7r,7r),r) at large imaginary times decreases 
albeit it shows no sign of singularity. This observation 



stands in agreement with results obtained at half-band 
filling^ and excludes the occurrence of a Kondo break- 
down. 

This evolution of the single-particle spectral function at 
J/t = 1 points to three distinct Fermi surface topolo- 
gies, sketched in Fig. pTJ In the paramagnetic phase 



(Fermi surface (1) in Fig. 



11| the Fermi surface consists 
(±7r, ±7r) points in the 



of hole pockets around the k 
Brillouin zone. Even though in our simulations charge 
fluctuations of the /-sites are completely prohibited, the 
Fermi surface volume accounts for both the conduction 
electrons and impurity spins. This Fermi surface topol- 
ogy maps onto that of the corresponding non-interacting 
periodic Anderson model with total particle density given 
by 1 -\- Uc and is coined large Fermi surface. In the 
antiferromagnetic metallic phase close to the magnetic 
order-disorder transition, (Fermi surface (2) in Fig. 11) 
the Fermi surface merely corresponds to a backfolding of 
the paramagnetic Fermi surface as expected in a generic 
spin-density wave transition. Here, a heavy quasi-particle 
with crystal momentum k can scatter off a magnon with 
momentum Q = (tt, tt) to produce a shadow feature at 
k + Q. It is only within the magnetically ordered phase 
that we observe the topology change of the Fermi sur- 
face to hole pockets centered around k = (±7r/2, ±7r/2). 
Due to the antiferromagnetic order and the accompanied 
reduced Brillouin zone, this Fermi surface topology sat- 
isfies the Luttinger sum rule. We coin it a small Fermi 
surface in the sense that it is linked to that obtained in a 
frozen /-spin mean-field calculation as presented in Eq. 
([9]). We have equally plotted the single-particle spectral 
function at higher values of J/t = 1.4 as a function of 
doping (See Fig. [T3|. For this coupling strength, and 
at half-band filling, the conduction band maximum is lo- 
cated at k = (tt, tt) with accompanying shadow bands 
at k = (0,0). The Fermi surface obtained upon dop- 
ing follows from a rigid band approximation, and yields 
the topology shown in Fig. 11 (2). The transition to 



the paramagnetic state shows up in the vanishing of the 
shadow features at k = (0,0). 

We note that our results are confirmed by simulations 
carried out using the larger cluster size Np = 4, see 
Fig. [M] In particular, we observe the small Fermi sur- 
face topology in the strong AF region with small coupling 
J/t and evidence for a large Fermi surface in the weakly 
ordered region {J/t > 1.2). 



C. Mean-field modeling of the single-particle 
spectral function 



We resort to a mean-field modelling in order to propose 
a scenario for the detailed nature of the topological tran- 
sition of the FS within the ordered phase. Within this 
model the topology change will correspond to two Lif- 
shitz transitions. Aspects of the single-particle spectral 
function can be well accounted for within the following 
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FIG. 12. (Color online) Single-particle spectra A^^(k,(jj) of the conduction electrons for J/t = 1 and t' /t = —0.3, at /3t = 80. 
The conduction electron density is reduced progressively from half-filling Uc — 1 (a) to Uc — 0.880 (h). The respective FS 
topologies are shown in Fig. [TT] 
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FIG. 14. (Color online) Single-particle spectrum A^^(k,(jj) of 
the conduction electrons obtained using cluster size Np = 4. 
confirming that the small Fermi surface topology in the lightly 
doped AF phase with coupling J/t = 1 is not an artifact of 
the smaller cluster results. The staggered magnetisation of 
the /-electrons is m{ = 0.657 =b 0.005. 
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FIG. 13. (Color online) Single-particle spectra A^^(k,uj) of 
the conduction electrons for J/t =1.4 and t' /t = —0.3, at 
/3t = 40. The shadow features present in (a) vanish as the 
occupation number ric is reduced. 



10 



mean-field HamiltoniarF^ 



H 



Ck+Qa 

ha 
/k+Qcr 




JV 
2 


A 



\ 



JV 
Jma 



A / 



(9) 



Here A, (/i) are Lagrange multipliers fixing the /- (c-) 
particle number to unity {ric) and m^, m( corresponds 
to the staggered magnetization of the conduction and /- 
electrons. The order parameter for Kondo screening is 
V (X {fJ^Ci^a- + c\ _^fi-(j). The large-A^ mean-field sad- 
dle point corresponds to the choice tti^ = ml = and 
y 7^ 0. As already mentioned this saddle point gives 
a good account of the hybridized bands we observe nu- 
merically in the paramagnetic phase. In the magnetic 
phase, one could speculate that Kondo screening breaks 
down such that the /-spins are frozen and do not par- 
ticipate in the Luttinger volume. This corresponds to 
the parameter set F = 0, m^ 7^ and m{ 7^ and 
leads to the dispersion relation of Eq. (Isl). Throughout 



the phase diagram of Fig. 11 the single-particle spectral 
function never shows features following this scenario. In 
fact, many aspects of the spectral function in the mag- 
netically ordered phase can be understood by choosing 
m^ 7^ 0, m( 7^ and V 7^ 0. This explicitly accounts for 
a heavy- fermion band in the ordered phase. It is in this 
sense that we claim the absence of Kondo breakdown 
within our model calculations. Fig. [15] plots the four- 
band energy dispersion relation £^n(k) obtained from the 
mean-field calculation in the magnetically ordered phase 
at t'/t = -0.3, J/t = 1 and at band fillings n^ = 0.898 
(Fig. [l5|a,b)) and Uc = 0.975 (Fig. [l5|c,d)) correspond- 
ing to the Fermi surface topologies labeled (2) and (3), 
respectively, in Fig. 
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In each case we have used the 
DC A value for the magnetization and varied V to obtain 
the best qualitative fit to the respective DCA spectrum. 
The change in the topology of the Fermi surface stems 
from a delicate interplay of the relative magnitudes of 
the magnetization rrig'^ and Kondo screening parameter 
V. 

In the DCA calculations, we are unable to study the 
precise nature of the topology change of the Fermi sur- 
face within the magnetically ordered phase since the en- 
ergy scales are too small. However we can draw on the 



mean- field model to gain some insight. Fig. 16 plots the 
dispersion relation along the k = (0,0) to k = (7r,7r) 
line in the Brillouin zone as a function of V. We would 
like to emphasize the following points, i) The transition 




(7r,7r) (0,7r) (0,0) 



(7r,7r) 




(0,0) (f,f) (7r,7r) (0,0) 

[kx, ky) 



(f,f) (7r,7r) 

[kx, ky) 



1 




10" 



10" 



10-3 

(7r,7r) (0,7r) (0,0) (tt, tt) 

[kxj rCy) 

FIG. 15. Mean-field band structures En(k) with J/t = 1 
{t'/t = -0.3). (a) V = 0.5, m^ = 0.066, m{ = 0.257, nc = 
0.898. The line intensity is proportional to the quasi-particle 
spectral weight, An(k,oj) . The Fermi surface topology (see 
closeup (b)) results from backfolding of the large- A^ mean- field 
bands. This corresponds to the Fermi surface (2) in Fig|ll[ 
(d) V = 0.260, m^ = 0.190, m{ = 0.693, Uc = 0.975. The 
Fermi surface topology (see closeup (c)) corresponds to that 
of hole pockets centered around k = (±7r/2, ±7r/2) as depicted 
by the Fermi surface (3) in Fig. 11 



from hole pockets around the wave vector k = (tt, tt) to 
k = (±7r/2, ±7r/2) corresponds to two Lifshitz transitions 
in which the hole pockets around k = (tt, tt) disappears 
and those around k = (±7r/2, ±7r/2) emerge. The Lut- 
tinger sum rule requires an intermediate phase with the 
presence of pockets both at k = (±7r/2, ±7r/2) and at 
k = (7r,7r). This situation is explicitly seen in Fig. 16 'b) 
at V = 0.314. ii) Starting from the antiferromagnetic 
state with hole pockets at k = (±7r/2, ±7r/2) there is 
an overall flatting of the band prior to the change in 
the Fermi surface topology. Hence on both sides of the 
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FIG. 16. (Color online) Topology change (a) of the Fermi 
surface as obtained from the mean- field modeling En(k.) of the 
DCA results {J/t = 1, t' /t = -0.3 m^ = 0.066, m{ = 0.257, 
ric = 0.898). As apparent on the scale of the closeup (b), 
hole pockets around k = (il=7r/2, =b7r/2) as well as around 
k = (tt, 7v) are present at V = 0.314. 



transition, an enhanced effective mass is expected. Ow- 
ing to point i) the effective mass does not diverge, iii) 
The two Lifshitz transition scenario suggested by this 
mean-field modeling, leads to a continuous change of the 
Hall coefficient. However, since a very small variation of 
the hybridization V suffices to change the Fermi surface 
topology one can expect, as a function of this control pa- 
rameter, a rapid variation of the Hall coefficient as it is 
observed in experiments^. 



ground state of the single particle Hamiltonian of Eq. [9l 
Both a variational quantum Monte Carlo calculation^^ ^' 
as well as a Gutzwiller approximation^ show a rich phase 
diagram which bears some similarities but also impor- 
tant differences with the present DCA calculation. At 
small doping away from half-filling, the magnetic transi- 
tion as a function of the coupling J/t is continuous and 
of SDW type as marked by the back-folding of the Fermi 
surface. Within the magnetically ordered phase a first 
order transition occurs between hole (AF/j,) and electron 
(AFg) like Fermi surfaces. At larger hole dopings, the 
magnetic transition is of first order and the Fermi sur- 
face abruptly changes from large to small (AFg). The 
calculations in Refs. I27H291 are carried out at t^ = 0. At 
finite values of t' = —0.3t one expects the AFg Fermi- 
surface topology to correspond to hole pockets centered 
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On the 



around k = (±7r/2, ±7r/2) as shown in Fig. 
other hand, the AF/j, Fermi surface arises from a back- 
folding of the large Fermi surface (Fig. fT3|. With this 
identification, our DCA results bear some similarity with 
the variational calculation in the sense that the same 
Fermi surface topologies are realized. However, in our 
calculations we see no sign of first order transitions and 
no direct transitions from the large paramagnetic Fermi- 
surface topology to the AFg phase. It is equally impor- 
tant to note that in both the variational and DCA ap- 
proaches no Kondo breakdown is apparent. In particular 
the variational parameter V which encodes hybridization 
between the /- and c-electrons never vanishes^^ s^. On 
the other hand, T = CDMFT calculations of the three- 
dimensional periodic Anderson model on a two site clus- 
ter and finite sized baths^'^ ^^ have put forward the idea 
that the magnetic order-disorder transition is driven by 
an orbital selective Mott transition or in other words a 
Kondo breakdown^. In particular even in the presence 
of spin symmetry broken baths allowing for antiferromag- 
netic ordering a big enhancement of the /-effective mass 
is apparent in the vicinity of the magnetic transition. In 
the paramagnetic phase, the low energy decoupling of 
the /- and c-electrons stem from the vanishing of the hy- 
bridization function at low energies. Our results, which 
makes no approximation of the number of bath degrees 
of freedom and which allow to treat larger cluster sizes 
do not support this point of view. In particular as shown 
in Ref. 8 no singularity in the quasiparticle weight is ap- 
parent across the transition. Equally, Fig. [T2]which plots 
the single particle spectral function at finite doping shows 
no singularity in Z(k = (7r,7r)) . 



CONCLUSION 



D. Discussion 

The evolution of the Fermi surface across the magnetic 
order-disorder transition has been investigated by other 
groups and especially in the framework of a Gutzwiller 
projected mean- field wave function corresponding to the 



We have presented detailed DCA+QMC results for 
the Kondo Lattice model on a square lattice. Our DCA 
approach allows for antiferromagnetic order such that 
spectral functions can be computed across magnetic 
transitions and in magnetically ordered phases. At 
the particle-hole symmetric point (half-band filling and 
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t' = 0) we can compare this approach to previous exact 
lattice QMC auxiUary field simulations^^ ^^. Unsur- 
prisingly, with the small cluster sizes considered in the 
DC A approach the critical value of J /t at which the 
magnetic order-disorder transition occurs is substantially 
overestimated. The important point however is that the 
DCA+QMC approximation gives an extremely good 
account of the single-particle spectral function both in 
the paramagnetic and antiferromagnetic phases. Hence, 
and as far as we can test against benchmark results, the 
combination of static magnetic ordering and dynamical 
Kondo screening, as realized in the DCA, provides a 
very good approximation of the underlying physics. 
As opposed to exact lattice QMC auxiliary field simula- 
tions, which fail away from the particle- hole symmetric 
point due to the so-called negative sign problem, the 
DCA+QMC approach allows us to take steps away from 
this symmetry either by doping or by introducing a finite 
value of t' . It is worthwhile noting that for cluster sizes 
up to 16 orbitals the limiting factor is the cubic scaling 
of the Hirsch-Fye algorithm rather than the negative 
sign problem which turns out to be very mild in the 
considered parameter range. 

Our major findings are summarized by the following 
points. 

i) We observe no Kondo breakdown throughout the 
phase diagram even deep in the antiferromagnetic 
ordered state where the staggered magnetization takes 
large values. This has the most dramatic effect at 
half-band filling away from the particle-hole symmetric 
point and at small values of J /t where magnetic order is 
robust. Here, magnetism alone will not account for the 
insulating state and the observed quasi-particle gap can 
only be interpreted in terms of Kondo screening. The 
absence of Kondo breakdown throughout the phase dia- 
gram is equally confirmed by the single-particle spectral 
function which always shows a feature reminiscent of the 
heavy- fermion band.»^ 



ii) The transition from the paramagnetic to the antifer- 
romagnetic state is continuous and is associated with the 
backfolding of the large heavy fermion Fermi surface, 
iii) Within the antiferromagnetic metallic phase there 
is a Fermi surface topology change from hole pockets 
centered around k = (tt, tt) at small values of the mag- 
netization ml to one centered around k = (7r/2, 7r/2) 
at larger values of ml. The latter Fermi surface is 
adiabatically connected to one where the /-spins are 
frozen and in which Kondo screening is completely 
absent .^^^^ This transition comes about by continuously 
deforming a heavy-fermion band such that the effective 
mass grows substantially on both sides of the transition. 
Within a mean-field modelling of the quantum Monte 
Carlo results the transition in the Fermi surface topology 
corresponds to two continuous Lifshitz transitions. As 
a function of decreasing mean-field Kondo screening 
parameter V and at constant hole doping, the volume 
of the hole pocket at k = (tt, tt) decreases at the 
expense of the increase of volume of the hole pocket at 
k = (7r/2,7r/2). Owing to the Luttinger sum rule, the 
total volume of the hole pockets remains constant. It 
is important to note that the energy scales required to 
resolve the nature of this transition are presently orders 
of magnitudes smaller than accessible within the DCA 
with a quantum Monte Carlo solver. 
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